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ABSTRACT 


This dissertation describes the design of an adaptive controller for single-input 
single-output (SISO) systems with guaranteed bounds on the transient response, and 
robustness with external disturbances and unmodeled dynamics. Developed from a 
current approach called “LI adaptive controller”, we show that by adding two properly 
designed low pass filters at the input and at the estimator we can control the transient 
response and the sensitivity of the overall system to external disturbances and 
unmodeled dynamics. Global stability of the overall adaptive system is mathematically 
proven under the assumption that the system is minimum phase (i.e., with the zeros of 
the transfer function in the stable region) and bounds of the systems parameters are 
known to the designer. 

The extension of this approach to non-minimum phase systems, such as 
systems with flexible appendages, is also considered. We show that a non-minimum 
phase plant augmented with a properly designed parallel system results in a minimum 
phase system. The augmenting system most easily comes from the inverse of a 
stabilizing Proportional-Integral-Derivative (PID) controller, designed to be least 
sensitive to parameter uncertainties. This approach is applied to a flexible arm in a 
testbed at the Naval Postgraduate School, called the Flexible Spacecraft Simulator 
(FSS), which emulates realistic conditions in space. Experimental results prove the 
effectiveness of the controller presented in this dissertation. 
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EXECUTIVE SUMMARY 


Adaptive control systems have been proposed to control systems with uncertain 
dynamics, widely varying under different operating conditions. They have been the 
subject of research for a number of years and stability has been established under ideal 
conditions. In spite of the amount of research performed on this class of control systems, 
sensitivity to external disturbances and modeling errors together with poor transient 
response performance made engineers skeptical to their applications. In particular, the 
Model Reference Adaptive Control (MRAC) approach, which seems to be the most 
effective under ideal conditions, can easily become unstable in the presence of 
unmodeled dynamics not accounted by the system model. Where instability means a 
complete loss of the platform, such as in space applications, adaptive control is far from 
being considered a possible solution to control problems. 

In this dissertation, we address some of these issues, such as transient response 
and robustness in the presence of external disturbances and modeling errors. The 
proposed adaptive control system is derived from a recently introduced adaptive control 
approach called LI Adaptive Control, with LI referring to the norm of a certain operator 
affecting the tracking error. We show that the addition of two simple low pass filters, 
one at the control input and one at the parameter estimator, eliminates high frequency 
components, which adversely affect both the transient response and the sensitivity to 
disturbances. These two additional filters must be properly designed in order to 
guarantee stability and fast convergence. 

A major drawback of this proposed scheme, relative to more traditional adaptive 
systems, is the need for a nominal model for the system and bounds on the parameter 
uncertainties. In other words, while transitional adaptive controllers assume the system is 
almost completely unknown (a black box), in the proposed scheme we assume the 
system to be within bounds (a gray box). In most, if not all applications, the gray box 
assumption is more realistic, since the order of magnitude of the dynamics is understood. 

A further issue examined is the extension to non-minimum phase systems. This 
problem is addressed by replacing the actual system with an “augmented” system, 
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implemented by adding a dynamie system in parallel. It is easy to show that a non¬ 
proper stabilizing eompensator, like a proportional-integral-derivative (PID) eontroller, 
ean be used by defining the inverse of its transfer funetion (thus proper) as the 
augmenting system. The eombination then beeomes minimum phase, and it ean be 
eontrolled by an adaptive eontroller. Applieation to a flexible arm shows the 
effeetiveness of this result, and testing at the Spaee Researeh and Design Center (SRDC) 
at the Naval Postgraduate Sehool emulates a flexible roboties arm freely moving in spaee. 

Although all stability results are proven mathematieally, a number of issues 
remain to be addressed. Speeifieally, further researeh is required to ehoose the optimal 
number of design parameters for best performanee of the proposed system. 
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I. 


INTRODUCTION 


A. MOTIVATION AND BACKGROUND 

The particular feature, which distinguishes adaptive control systems from other 
control techniques, is the adaptive loop that provides real time estimates of the Plant 
Dynamics. 

A general structure of an adaptive controller is shown in Figure 1. The system we 
want to control, called the plant, is modeled in terms of a set of parameters 9 *, unknown or 
partially known. The goal of the controller is to estimate the plant’s parameters and adjust 
the controller’s parameters accordingly. 



el 




Figure 1.1 Model Reference Adaptive Control. 
After [1] 
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As we can see in Figure 1, there are two loops. One is the standard feedback loop 
from the output j back to the input u . The other is the loop from the input-output of the 
plant to the compensator parameters. These two loops can interact with each other and 
cause instability in the system. 

This approach has been around for a number of decades, and stability under ideal 
conditions is well understood [1]. 

In spite of the seemingly attractive feature of self-adjusting to time varying 
dynamics, researchers and engineers have always been skeptical of adaptive control and, as 
of today, it is not widely used. There are a number of reasons for this skepticism. In 
particular: 

1. Because of internal modulation effects, fast adaptation produces a high 
frequency input. This input excites unmodeled dynamics easily, and causes 
systematic instability. Typical examples, from Rohrs [2], have shown that 
under very mild conditions of unmodeled dynamics that adaptive control 
system can easily become unstable. 

2. There is no guarantee on the bounds of the transient response of both input and 
output signals. Literature review [2-80] shows that most results do not provide 
any guidelines on transient performance analysis or explain how to improve 
transient performance. 

During the past few years, there has been greater interest in adaptive control, as 
evidenced by its prevalence in the literature. A particular approach by Naira 
Hovakimyan [1] shows that the addition of a filtering term in the loop greatly improves 
both robustness and transient response behavior. This approach is called “LI adaptive 
controller”, since it guarantees an arbitrarily small bound on the transient response 
using the small gain theorem on the L, norms of the error operators [81]. It was 
developed for state feedback and then extended to input-output systems. However, this 
extension does not follow the standard pole placement problem in classical adaptive 
approaches [82], and therefore, yields in a more general parameterization. This results 
in an adaptive system for a restrictive class of linear systems, but with the capability of 
compensating for a wider class of nonlinearities. 

The results shown in [81, 83-88] have proved the effectiveness of this approach 
in terms of robustness and transient behavior. 
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In this dissertation, we address the problems of designing a elass of Adaptive 
Control systems whieh yield fast adaptation, thus good transient response, and 
robustness to modeling errors. 

In addition, we show that this approaeh ean be extended to non-minimum phase 
systems, where the eurrent approaehes are unable to stabilize the system. Although this 
approaeh requires more knowledge about the system in order to eontrol it, it is still 
attraetive in oases where preoise eontrol is needed in the presenoe of model 
unoertainties. The typioal ease we analyze is the preoision eontrol of a flexible arm in 
spaoe. 

B. RESEARCH GOALS 

After introduoing the properties of a general adaptive eontrol system and the 
Model Referenoe Adaptive Control (MRAC) system, we examine the drawbaoks of 
these systems and use LI adaptive eontrol to address them. With the proposed LI 
adaptive oontroller, not only is the traoking error driven to zero, but the transient 
response is eontrolled, due to adaptation. The final result is a eontroller whieh, under 
the assumption of low disturbanees, yields fast adaptation to parameter variations. 

Another goal is to strengthen the robust ability of the adaptive system. After 
proper ehoiee of low pass filter, referenee model, and other important elements, the 
proposed adaptive eontroller demonstrates an improvement in the system’s ability to 
deal with the problems of modeling errors, nonlinearities and external bounded 
disturbanees. 

C. DISSERTATION ORGANIZATION 

This dissertation is organized as follows: Chapter II introduees the fundamentals of 
adaptive eontrol systems, speeifieally the direet Model Referenee Adaptive Control 
(MRAC) approaeh. Rather than repeating the detailed mathematieal derivations, whieh 
are found in the standard referenees, a simple first order example is presented, whieh has 
all the “fiavor“ of a full theoretieal presentation, but does not inelude the “fat”. Also in 
this ehapter, analytieal arguments and eomputer simulations are used to demonstrate that 
fast adaptation results in high frequeney signals and a possible loss of stability. 
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Additionally, using this same example, it is shown that the addition of a simple low pass 
filter in the loops guarantees good transient performanee and better robustness to 
modeling errors. 

Chapter III diseusses the general extension to a system of arbitrary degree. 
Typieally, the plant in a MRAC system must be minimum phase. This means that the 
plant to be eontrolled ean be unstable, and that the zeros of its transfer funetions must be 
in the left half side s-plane (in the “stability region”). Additionally, the use of a priori 
knowledge of plant nominal dynamies and its uneertainties is diseussed in the design of 
an adaptive eontroller, which is globally stable with a transient response having an 
arbitrarily small bound. This response is obtained provided the effect of external 
disturbances is negligible. 

Chapter IV addresses the minimum phase assumption, where the adaptive system 
is extended to a class of non-minimum phase systems. It is shown that a non-minimum 
phase system can be “augmented” to a minimum phase transfer function by adding a 
properly designed parallel system whose response is zero at steady state. Further, it is 
shown that if the system is stabilizable by a standard proportional—integral—derivative 
(PID) controller, then the transfer function of the added system is just the inverse of the 
PID controller. This is particularly attractive, since a properly design PID controller is 
generally robust in the presence of model uncertainties. However, such a system has the 
drawback that its transfer function is not proper, due to the derivative action of the 
controller. But, since its inverse is implemented in conjunction with the adaptive 
controller, its transfer function becomes strictly proper, and thus easily realizable. This 
chapter also discusses an experimental implementation of this approach on a flexible arm 
in the satellite Research and Design Center (SRDC) at the Naval Postgraduate School 
(NPS). In this experiment, vibrations due to highly damped flexible modes in the system 
are properly compensated for so that the position of the end effector quickly settles to the 
desired position. 
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Finally, conclusions and further research are presented in the final ehapter of this 
dissertation. Additionally, in order to better illustrate the proposed adaptive algorithm, a 
number of computer simulations and the results of applieations are shown alongside the 
theoretieal presentation. 


D. NOTATION 

In this dissertation, we define “ p ” the differential operator as 

px(t) = x(t) 


Also, the notation 


yit) = 


Bjp) 

A(p) 


u(t) 


with B(p ), A(p) polynomials in p , indieate the differential equation 

p^y(t) + a^p^^^ y(t) + ••• + % y(t) = biP^^^u(t) + • • • + b^u(t) 
with , bj the eoeffieients of A{p ), B{p ), respeetively. 

In all situations where the Laplaee Transform exists, the differential operator p can be 
also interpreted as the Laplaee variable 5. 
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II. INTRODUCTION TO ADAPTIVE CONTROL 


A. PREFACE 

According to Webster’s dietionary, to adapt means “to ehange (oneself), so that 
one’s behavior will eonform to new or ehanged eireumstanees.” This is exaetly the goal 
of an adaptive eontrol system, to self adjust to ehanging dynamies and environmental 
eonditions. A typieal example is the design of a high-performanee jet fighter or a fire 
helieopter [1], They both operate over a wide range of speeds, altitudes, and weights 
(before and after loading water for the ease of a lire helieopter) and their dynamies are 
nonlinear and time varying as deseribed by 

x{t) = A{t)x{t) + B{t)u{t) ^, 

y(t) = C(t)x(t) +D(t)u(t) 

where x(0) = Xo is assumed to be the initial state, and A{t),B{t), C{t), and D{t) are 
funetions of the operating time t. 

As the system goes through different operating eonditions, the dynamie model 
A{t) , B{t) , C{t) , and D{t) ehanges eonsiderably. The adaptive eontroller strueture 
eonsists of a feedbaek loop and a eontroller with adjustable gains as shown in Figure 2.1. 



Figure 2.1. Controller Strueture with Adjustable Controller Gains. 

From [1] 
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In this chapter, two different approaehes to adaptive eontrol systems are briefly 
introdueed. They are the Indireet Model Referenee Adaptive Control (Indireet MRAC) 
and the Direet Model Referenee Adaptive Control (Direet MRAC) approaehes. This is 
followed by a simple first order example to show the teehniques to analyze eonvergenee 
and stability. Finally, the approaeh presented in this researeh is introdueed for the first 
order example to demonstrate the main features of the proposed algorithm. 

B. MRC, INDIRECT MRAC, AND DIRECT MRAC 

MRAC is derived from the Model Referenee Control (MRC) approaeh. A linear 
time invariant (LTI), SISO plant of MRC is shown in Figure 2.2. 



Figure 2.2. Model Referenee Control. 

From [1] 

In MRC, full knowledge of the plant dynamies is assumed. The eontroller must 
be designed to meet the designer’s requirements, identified in the referenee model, as the 
desired Input-Output (I/O) properties of the elosed loop. The goal of the MRC approaeh 
is to find a proper feedbaek eontrol law so that the elosed loop matehes a referenee model 
with transfer funetion Wm(p). 

The feedbaek eontroller c{^6*^ is designed so that all signals are bounded and 
the elosed loop transfer funetion from r{t) to y{t) is equal to W^{p). This implies that 
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the object of the desired controller c{0*^ is to drive the tracking error 
^i(0 = 3^(0 “ Jm(0 close to zero. The requirements of the MRC are as follows: 

1. The zeros of the plant must be minimum phase. 

2. A good understanding of the plant and performance is required. 

If the zeros of the plant are non-minimum phase, the MRC may have an unstable pole 
zero cancellation thus causing instability. When the plant parameters are not totally 
known, an adaptive system based on estimates of the plant parameters must be developed. 
The estimated parameters can be obtained by direct and indirect approaches. The 
resulting control schemes are known as MRAC. A block diagram of a direct MRAC 
system is shown in Figure 2.3, and an indirect MRAC system is shown in Figure 2.4. 

The difference between a direct MRAC and indirect MRAC is that a direct 
MRAC updates the estimated controller parameter vector directly from input and output 
data. In indirect adaptive control, the controller parameters are computed on the basis of 
the estimated plant dynamics. In this dissertation, the approach is based on the direct 
MRAC. Accordingly, the concept of the direct MRAC is discussed in the next section. 
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Figure 2.3. Indirect MRAC. 

After [1] 


Input 

Comma] 

r(0 



Figure 2.4. Direct MRAC. 
From [I] 
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c. 


MODEL REFERENCE ADAPTIVE CONTROL 


1. Introduction of Model Reference Adaptive Control 

The general idea behind MRAC is to ereate a elosed loop eontroller with 
estimated parameters that ean be updated direetly from input/output observations. A 
standard, simple direet MRAC is shown in Figure 2.5. 



As shown in Figure 2.5, the output of the system y(t) is compared to a desired 
response y„(0 from a reference model. The control parameters are updated based on 
the error e(t) between y{t) and y^(t) . The goal of the MRAC is to update the 
estimated parameters 0{t) so that the plant response matches the response of the 
reference model VF,„(p) . In what follows, a first order example of an adaptive control 
system is presented and the methodology for determining stability and convergence is 
shown. 
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2. Theory of Model Reference Adaptive Control 

A simple example of a first order ease is used to introduce the theory of MRAC. 
Assume a plant with transfer function 

= ( 2 . 1 ) 

p-a 

where a is unknown, and it is not necessarily positive or negative. This implies that no 
assumption regarding the original plant’s stability is made. Equation (2.2) depicts the 
state space model as 


x(t) = ax(t) + uit) 
y{t) = x{t). 

Let the reference model have dynamics 

mx/ m z x/ /r\ 

where > 0 for stability. The reference output signal (t) represents the desired 
output the plant has to track. The goal of the controller in MRAC is to drive the tracking 
error e{t) = y{t)-y^{t) close to zero. In order to do so, the plant dynamics represented 
by Equation (2.2) are expressed in terms of the reference model as 


x{t) = -a^^x{t) + (a^ + a)x{t) + u{t). 
Eurther, the unknown parameter is defined as 


(2.4) 


0 = a^+a (2.5) 

so that Equation (2.4) can be written as 

x{t) = -a^x{t) + dx{t) + u{t). (2.6) 

This yields the structure of the control input, in terms of the model parameter 0 , as 
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u{t) = -6x{t) + v{t) 


(2.7) 


which yields the desired closed dynamies 

Kt) = -a„x{t) + v{t). ( 2 . 8 ) 

If we know the parameter a in equation (2.2), then we can calculate the control 
parameter in equation (2.5), and use the control input signal as in equation (2.7) to obtain 
the desired eontrol behavior. However, since we do not know the system dynamics, we 
use a time varying estimated parameter 6{t) is used to replaee the true eontrol parameter 
6. In this way, the eontrol becomes 

u(t) = -6{t)x{t) + v{f). (2.9) 

Now, the problem is how to compute the estimated parameter 6{t) based on input/output 
data. In order to do this, the eontrol input expression of Equation (2.9) is substituted into 
the plant model Equation (2.6) resulting in. 

x{t) = -a^x{t) + 6{t)x{t) + v(0 (2.10) 

where the parameter error 6{t) = 6- 6{t) is the difference between actual value 6 and 
estimated parameter 6{t) . Now if the tracking error e{t) is defined as 

e{t) = y{t) - yjt) = x{t) - xjt) (2.11) 

and Equations (2.10) and (2.3) are combined to relate the parameter error 0{t) with the 
tracking error e(t) , the result is 

e(t) = -a^e{t) + 0(t)x(t). (2.12) 

In order to provide a viable estimate, the parameter 0{t) is updated as follows 
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(2.13) 


6{t) = -6{t) = Yx{t)e{t) 

where / > 0 is an arbitrary eonstant. The basis for this is shown as follows. 

Define the Lyapunov funetion 

V(e{,),0(,)) = ^\e(t)f +iy-‘ |0(of (2.14) 

Differentiating V along the trajeetories of the systems (2.12), (2.13), and (2.14), we 
obtain 

y {e{t), e{t)) = e{t)i^-a^e[t) + 6{t)x{t)) + r'0{t){-rx{t)e{t)) (2.15) 

whieh yields 

V{eit),0{t)) = -aM^)f <0 (2.16) 

By standard argument this shows that the estimated parameter is bounded and the 
error goes to zero [1]. 

In this example, we ean also verify that an upper bound of the error e(t) ean be 
easily determined. In faet, sinee V(e(t), 6{t)) is non inereasing 

v{e{o),m)>y{e{t),m) ( 2 . 17 ) 

we ean write from (2.14) that 

i|e(0)|" + iy-' |#(0)f > \e{,f + iy-' |^(of. (2.18) 

By ehoosing zero initial eondition x^{0) = x(0) , i.e., e(0) = 0 we obtain 

i y-' |0(O)f > t|<-(0f +^r' |e(of ^ ^k(0f (2.19) 

for all t > 0 . This shows an upper bound for the traeking error as 

, ,2 ^(0)' 

\e{t)\ < - (2.20) 

y 
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for all t>0. From Equation (2.20), the relationship between tracking error e(t) , 


adaptive gain /, and initial parameter error 6’(0) is found. By increasing the gain /, the 
maximum error can be made arbitrarily small, thus achieving any desired tracking. 
Figures 2.6 and 2.7 show simulations for different values of the adaptation gain y. 


input u(t) 

2 p < ^ I I I I I r 


0 - 


-2 


-4_I_1_1_^_I_I_1_I 

0 5 10 15 20 25 30 35 40 


sec 

plant x(t) and model x^(t) outputs 



Figure 2.6. Response of Modeled MRAC Adaptive System with y = 50 . 
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/ 


input u(t) 


-5 L\_/f_ I _ I _ I_I _ ^ _ I _I 

0 - 5 10 15 20 25 30 35 40 

sec 



Figure 2.7. Response of Modeled MRAC Adaptive System with y = 1000. 


These simulations show that, by inereasing the adaptation gain y , the maximum 
traeking error deereases. Flowever, this eomes with a price. As shown in the simulation, 
the input signal exhibits a higher frequency component. In the next section, these high 
frequency components will be discussed, including their excitation of unmodeled 
dynamics, which cause the system to lose tracking or to become unstable. This is seen in 
Figures 2.8 and 2.9, where the addition of unmodeled poles at p = -2Q + j*5 and 

425 

p = -20 - 7*5 , with transfer function —r- , causes the system to either 

p'+40/7 + 425 

becomes unstable or loose tracking (Figure 2.9). 
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Figure 2.8. Response of Modeled MRAC Adaptive System with y = 50 . 
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Figure 2.9. Response of MRAC with Unmodeled Dynamics and y = 50 . 
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Figure 2.9 shows that, by increasing the adaptation gain y in the MRAC, the 
maximum tracking error may not decrease due to the unmodeled dynamics thus leading 
to loss of tracking. 

3. Adaptive Control with a Filtering Action 

The arguments presented in Chapter II. C.2 show that there between faster 
adaptation (large y) comes at the expenses of a high frequency transient response 
which can cause instabilities in the presence of modeling errors. This can be seen 
from the input signal, repeated here for convenience as, 

u{t) = -0{t)x{t) + v{t). 

Using the definition of parameter error 6{t) = 6-6{t) and its relation to the tracking 
error e{t) shown in Equation (2.12), the input signal can be expressed as 

u{t) = -Ox{t) + v{t) +e{t)-\-a^e{t). (2.21) 

While the error e{t) can be made arbitrarily small by increasing the adaptation gain y , 
its derivative e{t) not only is not necessarily small, but it can become quite large due to 
the high frequencies in the transient response. 

A better overall response is obtained by introducing a Low Pass Filter with 
transfer function Cp {s) at the input signal, which yields 

u{t) = Cp{p)[-e{t)x{t) + v{t)) (2.22) 

This filter’s bandwidth should be larger than the reference model’s, so that the filter stops 
the high frequencies generated by the adaptive system without affecting the dynamics in 
steady state. 

The input signal u{t) can be written again in terms of 6{t) and e{t) as 

u{t) ^ Cp{p){-ex{t) + v(0) + Cp(p){e{t) + a^e(t)). (2.23) 


18 



The rightmost term shows the effect of the adaptation error, and it can be written as 

u{t) = [p + a)Cp{p)e{t) (2.24) 

where the transfer function (5 + a)C^(5) is proper (i.e., no differentiation) and stable. 
This implies that small values for e{t) also result in small values for u{t) , and the 
controller is always close to an ideal controller with parameter 6. 

Figures 2.10 and 2.11 show this behavior. The system is the same as the system 
presented in Chapter II. C.2, but it is observed that an increased adaptation gain y yields 
faster response, with better transient behavior. This adaptive control with a filtering 
action approach has been called LI adaptive control [81, 83-88]. 


input u(t) 



0 5 10 15 20 25 30 35 40 

sec 

plant x(t) and model x^(t) outputs 

-r-1-i-!-1-1-1-1- 

-plant x(t) output 

» - -model Xm{t) output - 



- 0.5 L- 1 -1-^-1-^-1-^-J 
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Figure 2.10. Response of LI Unmodeled MRAC Adaptive System with y = 50. 
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Figure 2.11. Response of LI Unmodeled MRAC Adaptive System with y = 1000 . 
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III. LI ADAPTIVE CONTROL FOR MINIMUM PHASE SYSTEMS 

A. OVERVIEW OF LI ADAPTIVE CONTROLLER 

In this chapter, the adaptive control formulation is extended to a general SISO 
system. In particular, it is shown that the plant can be parameterized in terms of the 
eontroller coeffieients, so that they can be estimated directly from the input and output data. 
Further, the effeetiveness of adding a LPF in the controller is addressed. By proper choice 
of this LPF, it is shown that the magnitude of the transient response due to adaptation can 
be made arbitrarily small for both input and output signals. The cost is a perturbation of the 
reference model, which is still globally stable. 

This chapter is organized as follows. Model definitions are introduced in 
Chapter III. B. Parameter estimation is given in Chapter III. C. Then, all properties and 
elements are combined to show the overall model in Chapter III. D. In Chapter III. E, 
the LI adaptive controller is analyzed under ideal conditions, without any disturbances 
or modeling error. Then, in Chapter III. F, analysis of the ideal case is extended to the 
eonditions of external disturbance and unmodeled dynamics. Finally, in Chapter III. G 
some examples are considered that compare MRAC systems and verify the theorem of 
LI adaptive controller. After this, a deeper exploration of the properties is possible. 

B. MODEL DEFINITIONS 


In the SISO cases, consider a system with dynamics described by 



y(t) = K m (0 + >^ o (0 

A{p) 

(3.1) 

where 

^ip) ~ P^ +a^p^ * +... + a^ 

(3.2) 


B{p) = p^ +b,p^-^+... + b^ 

(3.3) 


are monic. Then the reference model is defined as 
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(3.4) 


where 


Jmit) 


1 

D{p) 


v{t) 


D{p)=p^-^ +dy-^-^+...+d 


N-M 


(3.5) 


is a Hurwitz polynomial with degree 

dD{p) = dA{p)-dB{p) (3.6) 

and v{t) is an arbitrary external bounded input. 

In the aforementioned diseussion, and in most of the related literature, the 
assumptions on the plant are as follows: 

a. The degree of the numerator B{p), M, and the denominator A{p), N, are 
known. 

b. The sign of the high frequeney coeffieient kg is known (assume ko > 0 
without loss of generality). 

e. The high frequeney eoeffieient kg and the parameters kj,..., 
of the transfer funetion of the plant are unknown. 

d. The plant ean be stable or unstable. However, the numerator B{p) is 
Hurwitz, i.e., it has roots with negative real parts. 

1. Plant Parameterization 

Now, the goal is to parameterize the system in terms of the eontroller parameters. 
In partieular, following the pole plaeement approaeh, and assuming the following state 
spaee model, as follows: 

x{t) = Ax{t) + Bkf,u{t) 

(3.7) 

y(t) = CpX{t) 
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Using an observer and state feedbaek, the elosed loop poles ean be arbitrarily 
plaeed as shown in 

x(0 = (Ap -K^C^)x(t) + B^kou(t) + K^y(t) 
koM(t) = -LpX{t) + v{t). 

If the eontroller gain is ehosen so that 

det(^/ - +5^L^) = B{p)D{p) (3.9) 


and the observer gain as 


dQi{sl-A^ + Kf^) = P{p) 


(3.10) 


with P{p) arbitrary, then the Hurwitz elosed loop dynamies beeomes 


y{t) = 


B{p) 

D{p)B{p) 


v(0 


1 

D{p) 


v{t). 


(3.11) 


In order to determine a suitable plant parameterization, the eontroller in Equation 
(3.8) is written as 

k,u(t) = -L^ (si-A^+ y B^u(t) -Lysl-A^+ y y{t) + v{t). 


This expression ean also be written as 


kpA{t) = + v{t) 

“ P{p) Pip) 


where 


' ’ ' I’ p(p) p(p) 


u P P P’ P 


(3.12) 


(3.13) 

(3.14) 
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Define the veetors of eoeffieients 


and the filtered signals as 


h [/Zj,..., /z^ ] , 




te] 


p{p) 


N-2 

P 

P 

1 

yN-l 

N-1 


P 

1 


u{t) 


yit) 


then the plant model ean be written as 


1 


Dip) - 


Sinee the bounds on L are assumed to be known as 




equation (3.19) ean be expressed as 


1 


J(0 = 777- + (^0 - ) W(0 + (0 + ^ V, (0) 

Dip) 

In summary, the following statement can be made: 

Given the plant 

Aip) 


(3.15) 

(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 
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and an arbitrary Hurwitz polynomial D{p) with degree dA(p) - dB(p ), the plant can be 
represented as 


Dip) 


(3.21) 


where 


(pit) = 


u(t), 


N-\ 


Pip) 


uit), 


1 


7V-1 


-u(t) 


Pip) Pip) 


yit),- 


Pip) 


yit) 


(3.22) 


and 6^ e ^ 2 w+i controller parameters. 

It can be easily seen from Equation (3.21) that, in order to track the model ——, 

Dip) 

the controller has to be of the form 

kjiit) = -^(l)it) + vit) i3.23) 


C. PARAMETER ESTIMATION 

In this section, the problem of estimating the parameters P from the plant 
dynamics 

yit) = - j— ikMO + 0^<Pit)) + wit) (3.24) 

Dip) 

is addressed. In order to do so, let C^ip) be a transfer function such that the Cj,ip)Dip) 
is proper and is of the form 


CEip) = 


(5 + COe) 


(3.25) 


with L> N -M . Thus equation (3.24) can be arranged as 

Ce ip)Dip)yit) = Q ip)ikjiit) + (ko - k^)uit) + it) + k^p^ it)) + it) (3.26) 


or 


CEip)Dip)yit) = C^ip)ikjiit) + ^ (pit)) + w^it) 

25 


(3.27) 



where Wj,{t) = Cj,{p)D{p)w{t). If Wj;{t) is set as the upper bound of w^(t), then the 
parameter of equation (3.26) and (3.27) ean be written as: 

CEiP)Dip)y{t) - k^c^{t)u{t) = e^X (0 + w^{t) (3.28) 

where 6^ e ^ 2 w+i 

Notiee that the parameter veetor 6 appears in (3.28). Therefore, in order to 
estimate it, eall ^(t) an estimate and define the error as 

e{t) = Ci,{p)D(p)y(t) - k„p^{t)u{t) - 0'^ {t)X{t) (3.29) 

This ean be easily related to the parameter error 6{t) = 6-6{t) 

e{t) = e'^{t)X{t) + w^{t). (3.30) 


Based on this equation, the parameter 6{t) ean be estimated as follows: 

z,(0 


0.(t) = p,(t) 


cr w^(t) + X(t) 


-e(t), i = 0,...,2N 


(3.31) 


with 


0 , if 6»,.(0 = e.^ and XXt)e{t) > 0 


Pi (0 = T 0, if 4 (0 = 6». ^ and X,. {t)e{t) < 0 


^ p , otherwise 


(3.32) 


and p > 0 arbitrary. By this definition, the estimated parameters 6^{t) is bounded as 
shown in Figure 3.1. When the estimated reaehes its maximum (minimum) value, 
adaptation is turned off (//;(t) = 0) as long as it tends to inerease (deerease). 
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Figure 3.1. The Figure of Estimated Parameter 0^ (t). 


From the aforementioned information, it can be shown that for all t>0 that 


r 


a- o.^^<ext)<e. 


< 


b. 


o'" {t)x (0 < (0 (cr IWg (0| + ||x (Oil) 


c. 


0{t) 


<jUS2it) 


V 


with (0, ^2 (0 > such that for all t, T > 0 


(3.33) 
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^ "m"" 


t+T ^ j 

j|^,(r)|jr< < 


2 // 


m 


if Wg (r) = 0 for all r 


(3.34) 


-t —otherwise 


2^ a 

The proof is shown in Appendix A. 

In the absence of disturbances, (t) = 0 for all t > 0 . Therefore, Equation 
(3.31) can be restated as 


= i = 0,...,2N 

1^(01 


(3.35) 


with 


^ 0, if 0,{t) = e.^Ax and X^{t)e{t)>Q 


A, (0 = ^ 0, if e. it) = A. ^ and A,. {t)e{t) < 0 


(3.36) 


jj ., otherwise 


In the case of no disturbance, in order to avoid singularities during 
implementation, the denominator ||A(t)||^ can be replaced with a + ||A(t)||^ for any 
arbitrary constant a , without loss of generality. In developing the simulated model, 

II l|2 

the arbitrary constant is set to one, a = 1. Thus, the denominator becomes 1 +1| A(t)| . 

D. OVERALL MODEL 

As in the first order system example presented in Chapter II, the filter {p) is a 
low pass filter with a DC value of one, C^(0) = 1. This filter is critical to the stability 

and performance of the overall system, and it is addressed in greater detail later in this 
section, including the criteria for its selection. 

Let the input of the adaptive control system be represented as 
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kji{t) = C^(p)(-^^(0^(0 +v(0) 


(3.37) 


where v{t) is the external input and is as previously described in Equation (3.22). 
The estimated parameters 6^{t) are computed, as shown in Equation (3.31). 

Now the control input in (3.37) is combined with the parameterized plant model 
in Equation (3.24) to obtain the closed loop dynamics model as 

yit) = ^^0' mt) + v(0)++ Mt) 

D{p) D{p) ^ ^ (3.38) 

kji{t) = Cp ip)[0^ it)(pit) + v(0) - Cpip)[0^(pit)) 

where 0^ = k^ ^ is the true parameters vector. Again 0{t) , the estimate, 

is computed as in Equation (3.31), and 0{t) = 6-6{t) is the parameter error. 



eg{t) 

Eigure 3.2. Dynamics Model of the Adaptive System. 


In order to better understand the global stability of the overall adaptive system, 
equation (3.38) can be represented by the dynamic model in Eigure 3.2 where we define 
the error term 
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eg{t) = Cp{p)(e^{t)^{t)) 


(3.39) 


as the effect of parameter error. The overall model (p) is obtained by combining 
equation (3.38) with the two plant parameterizations 

A{p) 

and 


Dip) 

where ^K<kM . This yields the dynamics model 


u[t) 


= {p)Cpip)(0^it)^it) + v(0) + H^{p)wo(t) 


(3.40) 


defined as follows: 

yit) = {p) yit) + H^^{p)Cpip)i^d^it)^it)+vit)) + H^^{p) w^it) 

kjiit) = [p) yit) + //„„ {p) Cpip)[e'^it)<f>it) + v(0) + H^^{p) Wait) 


where 
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H,,(P)=(1-C,(P)) 1- 


A(p) 


B{p)D{p) 


\+yC^[p) 




(1 + /) 


[\ + yCp[p))D{p) 




(l-C^(p)) f A{p) ^ k{p) 

(l + rCp{p))[^B{p)D{p) P{p)D{p) 


\ 

1 +1 

y 


H^P) 

hM 


1 

\+rc^{p) 

CF{p)D(p) f A{p) ^ 

\ + yCp{p) yB{p)D{p) ^ 

Cp{p)D{p) ( A{p) ^ kjp) 
l + [.Bip)D(p) P(p)D(p) ^ 


k{p) 

P{p)D{p) 


(3.42) 

k k 

and y = — -\ in the interval 0<y< ^ -1. The Proof is shown in Appendix B. 

k^ k^ 


The system H^^^{p) in (3.40)-(3.42) represents a perturbation of the dynamie 
model (3.24) due to the presenee of the Low Pass Filter Cp{p). In faet, if we let 
CpiP) P ’ we assume a filter with infinite bandwidth the terms in (3.42) 

yield ( negleet the disturbanee terms for eonvenienee) 


»hp)=(i-c,(AT 


A{p) 


B(p)D(p) 


\ + yC^{p) 


and 


Hyo{p) = 


(1 + y) 


{\ + yCp{p))Dip) D(p) 


This yields 


y(t) = tJ— (0^(0+v(0) 

Dip) 


(3.43) 


(3.44) 


(3.45) 


31 



which is the standard dynamic operation in the presence of parameter error 6{t) . For the 
input term, the analogy is slightly more complicated, but with Cp{p) = 1, simple algebra 
yields 

Ku{t) = -e\t)m + v{t). (3.46) 

Now the issue is how to choose the filter Cp{s) so that the transformation 
Hr^f{s) is exponentially stable. Examination of Equation (3.42) yields two issues: 

1. (l + yCp ( jco)^ * must be stable for all values of / in the interval 

^/max ■ 

2. must be stable for all allowable /, D{ja)), 5 ( 7 ®), and 
A{jo)). 

This leads to the following two conditions. 

Condition 1. Eet (p) be of the form 


Cp{p) = 


r \ 

f L 3 

V E + y 

[{p + CD,fj 


Coip)Ciip) 


(3.47) 


where L> N -M the relative degree of the plant. Also, let o, . Erom the Bode plot, 
shown in Eigure 3.3, it is always possible to choose (U, and with gain margin larger 
than the parameter uncertainty 7 . 
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Cp{co)\ 


► CO 



Figure 3.3. Bode Plot of Filter Cp{p). 


For the stability of (^7 - 77^^ (p) j , recall the expression for 77(p) 


H„(p) = (l-C,(p)) 1- 


Mp) 


B{p)D{p) 


\ + YCp{p) 


Since A(p) , B{p) , and 7)(p) are all monic polynomials and B{p)D{p) is 
Hurwitz by assumption, it can be seen that 


bm 

p^cc 


\ Ajp) ^ 

V B{p)D{p)j 


= 0 


(3.48) 


Since A(p) and [B(p)D{p)) have the same degree, this implies that 


\ A(p) ^ 

V B{p)D{p)j 


IS 


low pass. At the same time, with ( 5 ) being a low pass fdter with unity DC gain, this 
implies that (l-C^(p)] is a high pass filter. Therefore, the term 
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(i-c,(p)) 1 


A{p) 


B(p)D(p) 




JJ 


can be made arbitrarily small by inereasing the 


bandwidth of C^i^p). This is summarized by the following eondition; 


Condition 2. For every positive eonstant we ean always find (Ug sueh that 


max< 


JO) 


1 -- 


AUo)) 


B{j(o)D{ja>) 


= <1 


(3.49) 


® [U® + ®o| 

for all CO. 

These two eonditions are at the basis for ehoosing the filter 

f V ^ ^ 

CAp)= (3.50) 

\P + G)a)y[p + Q}^) j 

given any arbitrary positive eonstant < 1, and ehoosing cOq to satisfy Condition 2, 
(i.e., Equation (3.49)). Then ehoosing o,, so that the gain margin is 


< |l VMAx\dB (3.51) 

where co^ is the phase erossover frequeney of the filter, guarantees that both Conditions 1 
and 2 are satisfied. 

E. ADAPTIVE LOOP FEEDBACK FUNCTION 

The parameter error term in Figure again defined as 

eg{t) = Cp{p)(^{t)(j^{t)) 

This relates the parameter error 6{t) with the error eg{t) and it represents a 
transformation Hg : (p{t) eg{t) . 
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Recall the key results concerning parameter estimation previously presented in 
Chapter III. C. In the case of no external disturbances, the adaptation algorithm yields 
the following upper bounds 

0^{t)X{t)\<\£,{t)\\\X{t) 
and 

0(t) <//^2(0 

where 

X(t) = C,(p)m (3.54) 

Both terms £^{t) and are in and their magnitudes are upper bounded as 

r, ^(0) ' 

j|f;(r)|jr< - (3.55) 

The goal of this is to show that all the bounds above (3.52)-(3.55) yields an upper 
of the error e^(t) due to adaptation as 



\egit)\<P{t)sup\\(piT)\\ 

0 <T<t 

(3.56) 

where 


lim /3(t) = 0 

r—>00 

(3.57) 

and 


fcoY 

\m\<K — 

(3.58) 


where is a constant depending on the relative degree L, Og is the bandwidth of the 
fdter Cp{p ), and cOj, is the bandwidth of the fdter Cj,{p). 

This is very important since, going back to Figure 3.2, equation (3.56) shows that 
the effect of adaptation eg{t) is bounded by a feedback gain /3{t) , which not only goes to 
zero (3.57), thus making the system globally asymptotically stable, but also its maximum 


(3.52) 

(3.53) 
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value (3.58) depends on the bandwidths of the two filters and it ean be made arbitrarily 
small. The rest of this seetion is devoted to prove this result and, not surprisingly, is very 
teehnieal in nature. 

The abounds (3.56)-(3.58) will be shown as sequenee of faets. 


Fact 1: 


ee(t) = 


kCe(p)j 


eAt) = 


kCe(p)j 


Co(p)(e^(t) + e^(t)) 





< 

V 


(3.59) 


0<A<t 


1^1 (0| + p\ 1 ^ 2 (^)| Sl(^-^)) dT sup ||^/)(;l)|| 

g ^ 0</l<r 

where ^j(t)and determined by Equations (3.33) and (3.34), and 


S, 


1 

= ‘■"‘'"'IS 


fL-l 


\k=o kl j 


(3.60) 


Proof : First, rearrange Equation (3.39) 

eg{t) = Cg(p)Ci(p)(^^(0^(0) ^ 
Then the error ean be written as 


^ C,(p) ^ 

\CAp)j 


c,ip)[c,ip)(0At)m) 


eg(t) = 


^C,{p)^ 


c, 


,(p){(^^(0C,(p)^(0) + [c^(e)(^"(0^(0)-(^"(0C,(;?)^(0) 


\Ce(p)j 

If the term in the eurly braeket is defined as e,{t) , this expression beeomes 


eeiA = 


^C.ip)'^ 


eAt)- 


\Ce(p)j 

The term c^(t)isthe sum of two terms 


(3.61) 
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The first term ean be written as 

(0 = Co (;?) ( O'" {t)Ci, {p)^{t)) = C^(p)e^ (t) (3.62) 

with 

(0 = e\t)c,{p)m = e\t)x{t). (3.63) 

Additionally, the second term can be expressed as 

eJO = Co(p)[c,(p)(^^(0^(0)-(^"(OQ(p)#))] = Co(p)e^(0 (3.64) 

with 

(0=c, {p)[e^ {t)m) -[o^ it)c, {p)m\ (3.65) 

From equation (3.33) we bound (3.63) as 

(0| = \e\t)C,{p)m\ = \0Ht)X{t)\ < |^,(0|||3^(0|| (3.66) 

Since X{t) = C^{p)(/){t) , the following expression can be written: 

dr sup 

J 0</l</ 

By the way we choose the fdter Cj,{p) in (3.25), the impulse response (t) > 0 
for all t, so the rightmost integral is C^(0) which is equal to one. Therefore, X{t) can be 
bounded as follows: 

|A(t)| < sup ||^zi(yl)||. (3.67) 

0<yl<f 

Combine this with (3.66) to obtain the upper bound 

^^(0 ^|^i(0|||3f(0||^K(0|sup||^(2.)||. (3.68) 

0<x<t 

This shows the bound of the Fact, in (3.59). 

For the bound on (t) , expressed in Equation (3.65), it is shown in Appendix C that 



1 ^( 0 | ^ j \ce {t - 2 ')| ||^(^)||^ j \ce (^) 

0 Vo 
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W j C^(T)^(^t-T^dT dX 


(3.69) 


The derivative ^(t) is upper bounded as in (3.53), so that (3.69) above beeomes 


eAt) <h\\^ 2 A)\ f c^{T)dT dX sup ||^(ro) 


(3.70) 


where again C£(t) > 0 is the impulse response of the fdter = 


{p + g>e) 


-—- . Then 


+00 L L-1 

J cAT)dT= j e-^^AT = g^(co^{t-T)) 

t — X t — X V / 


(3.71) 


This shows the seeond inequality in (3.59), and Faet 1 is proved. 


Fact 2 ; The error signal 


QED 


ycEip)) 


(3.72) 


ean be upper bounded as 


Xt) <A0sup||^(r)|| 


(3.73) 


where 


lim y0(t) = 0 

t^oo 


p \ COq 


m < ^+- ^ p(0) ,forallt>0 

I 2)1// I" II 


(3.74) 
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4 = 


+“lL-l 

A 




Proof : From Equation (3.59) 


-2x 1 

e ax 


r I.r - 


1 {k,+k^)\ 


zz 




eg{t) = 


\Ce(p)j 


eAt) 


with 


Using the bounds in Equation (3.59) yields 


(0 ^ A(0sup||^(A)|| 


0<A<t 


Where, using Schwartz inequality, 


lAWhj 


= I 


-a>„{t-T) I 


^2 f +® 2 ^2 


?j(r)Jr<Oo j e J|fj(r) dr 


This implies that 


And, therefore 


lA ( 0 |s 




f 2 

®0 ^-2a,„T 


2 f 

^( 0 ) ' 

^-2co, 

0 7 

V 

2m j 


for all t > 0 . Eurthermore, sinee (t) e , then 

lim y0j (t) = 0 

r->-oo 

Therefore, Equation (3.75) ean be rewritten as 



(3.75) 


(3.76) 


(3.77) 


(3.78) 


(3.79) 
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For the second part of eg(t), recall from Equation (3.59) 


with 


(0 ^ A(0sup||^(r)|| 


0 <r<t 




(3.80) 

(3.81) 


Since s^it) e , then 


= 0 . 

/->OD 


Again using Schwartz inequality in the convolution integral in (3.81), we obtain 
an upper bound for (t) 


2 f +-X 




y V 0 


< A,. 

L \ 


I A 

2a>r 


m 


(3.82) 


Since Cg(t) = Oge , for t > 0, it is easy to see that 


A(0 = Co(0*A(0 


(3.83) 


limy02(O = O 


still yields the limit 
and the same upper bound 

1 

/ +00 ^2 f' 


<A, 


A 


\2o)e 


^(0) 


(3.84) 


(3.85) 


Combining Equations (3.78) and (3.85) yields the result in Eact 2. 

QED 

As noted, the dimension of the adaptive gain, // > 0 , is inverse time or hme ', as 
demonstrated in Equation (3.31). As a consequence, — in Equation (3.74) is 

Q) 

dimensionless. 

The bound on the term e^{t) can be found in Equations (3.73) and (3.74) as 


k/0 ^ AO sup|Ar)| 

0<r<r 


(3.73) 
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, for all t>0 


(3.74) 


m< 




m 


From Equations (3.73) and (3.74), e^(t) can be minimized by ehoosing the adaptive gain 
ju in Equation (3.31) in terms of the parameters co^, co^, and degree L of the fdter 
Cp ( 5 ) . And then, by straightforward differentiation the smallest upper bound is obtained 


// 


whieh yields the bound on the adaptive transformation Hg(p) 


eAt) <K 




V J 


sup (/){t) 


0 <r<r 


(3.86) 


(3.87) 


for some eonstant kp , dependent only on the filter’s order and degree, L> N-M . 


This result shows how the upper bound on the transient response e^(t)due to adaptation 
depends in a fairly simple fashion on the ratio of the bandwidths and cOp of the two 


filters Cp{s) and Cp{s) in the eontrol and the estimation loops. As 






0 the effeet 


of adaptation in both input and output signals tends to zero 

In the next seetion, we show this result formally, based on the what has been 
shown in this seetion, summarized below: 




C,{p) 

Ce{p) 


eA^) 


it ean be shown that eg{t) has the same bound as e^{t) 


(3.88) 
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f 


< 


V 


AOsup||^(r) 

0<r</ 


m<K 


^ CO. ^ 


]_ 

4 




lim P{t) = 0 

^->00 


(3.89) 


F. STABILITY AND BOUNDEDNESS OF THE ADAPTIVE SYSTEM 

In this section, the concepts up to this point are combined and culminate in the 
main result of this dissertation. 

First define the vector z{t) 


zit) = 


yit) 

u{t) 


(3.90) 


Assuming zero disturbances, from (3.40) and (3.89) we obtain 


r z{t) = H^^^[p)[eg{t) + Cp{p)v{t)) 

i (3.91' 

[ eg{t) = Hg{p)(^^{t)) 

where H^^j(^p) is exponentially stable, and Hg{p) is bounded by /3{t) in Equation 
(3.89). Also, recall that (/){t) is a vector of filtered input and output signals, which can be 
written as 

m = ^{p)m (3.92) 

where ^(p) is exponentially stable transformation in (3.17) and (3.18). 

The whole system, obtained by combining Equations (3.91) and (3.92), is shown 
in Eigure 3.4. 
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Figure 3.4. Dynamics Model of the Adaptive System with No Disturbance. 

Since, ^ 0 (Equation (3.89)), and and 0(p) are exponentially 

stable, the closed loop system is exponentially stable and all signals are bounded. 

Also, if the reference input and output signal are defined as 

Zrefit) = H^^f{p)Cp{p)v{t) (3.93) 
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From these, it is easy to see that 


^IL ^ \^refip)\^ |®(i^)L (II^IL +lk-/|L) (3.99) 


whieh yields 




(3.100) 


with 


Inequality (3.100) is valid only when the gain K is smaller than one as 0 < .ST < 1. 
From the results above, reeall 


(3.101) 


Hg{p)^ = sup /3{t) <k. 


0<t<co 




V y 


(3.102) 




where depends on the filter order only. Therefore, it is observed that as —^ q , the 


0)^ 


overall loop gain K goes to zero and, therefore, the error terms 


(3.103) 


are sueh that 


and 


lim sup y{t) - y _ (0=0 (3.104) 

0</<+a) 


lim sup w(0-Mre/(0 =0. (3.105) 

^—^00 0 <f<+oo 
COe 


In other words, the maximum value of the transient response due to the adaptive 
eontroller tends to beeome arbitrary small. 
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In the next seetion, we show that eomputer simulations support the results 
presented above. 

G. EXAMPLES OF APPLICATION 

I. LI Adaptive Control in Ideal Case (No Disturbance) 

Consider an unstable plant 

yit) = k,—u{t) (3.106) 

p-l 

where 1 <k,< 4, and the referenee model 


p + 2 


Then 


B(p) = l ,A(p) = p-l, D(p) = p + 2. 


(3.107) 


(3.108) 


First, the true coefficients are determined, assuming k^ =3 as the true value of the 
parameter. From A(p) and B(p ), the related space state equation can be determined and 
is expressed as follow: 

x(t) = A^x(t) + B^u(t) ^ 

y(t) = CpX(t) 


where 



Following the steps mentioned above, the closed loop system is obtained 

-^(0 = -S^L^)x(t) + B^v(t) 

y(t) = C^x(t) 


where controller gain =3 is derived from Equation (3.9). 


45 



Let the observer be 

x{t) = A^x{t) + B^u{t) + K^{y{t) - y(t)) 
y(t) = CpX{t) 

(3.111) 

where P{s) = p+ 2 and the observer gain = 1 is derived from Equation (3.10). 

From Equations (3.15), (3.16), (3.22), and (3.23), the system ean be represented 
as 

kQu{t) = -^[kQ-k^) k^~^^{t) + v{t) = -0^^{t) + v{t) (3.112) 


where 


0^ =[2 3 4.5f 

(p{t)= u{t) - -u{t) - -y{t) 

^ s + 2 s -\-2 


(3.113) 


Now, with ^0 = 3, a suffieiently large value for Og must be ehosen. Additionally, 
aeeording to Condition 1 and Equation (3.49), the magnitude of the response of 
^ A(p) ^ 

1-must be less than 1. The magnitude of the Bode Plot of 

I B{p)D{p) 


\ Mp) ^ 

V B{p)D{p)^ 


is shown in Figure 3.5 and it is smaller than one (negative dB’s) for all 


frequeneies larger than than 2.5 rad/see. This ean be used as a guideline to determine the 
bandwidth co^ of the fdter Cp{s). 
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bode of (1-(A{s)/(B(s)D(s))) 



Figure 3.5. The Bode Plot of 


A{p) ^ 

V B{p)D(p)j 


With a choice of Oq = 5, the LPF becomes (p) = 
must verify that Condition 2 is satisfied. Recall this condition 


yP + a„j 


s + 5 


and we 


max max \h ( jco)] < 1 

CO o<r<rMAx' ' 


where 




A{p) 


B{p)D{p) 


(p)y 

For a number of different realizations of y, is plotted, as shown in 

Figure 3.6. The simulation results conform to the requirement of Condition 2. 
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bode of Hyy(s) 



Figure 3.6. The Bode Plot of for 100 Plant Realizations. 


Next, if cOj; is chosen such that cOj^xo^, and the degree L of Cj;{p) is guaranteed to 
meet L> N -M , then Cj;{p)D{p) will be proper. In this example, let =100 rad/sec, 
L = l = M -N, and, therefore. 


Ce{p) = 


co^ 


100 


(P + ®£) (p + 100) 


(3.114) 


Since this case is implemented with no disturbances, W£(t) = 0 . Figure 3.7 and Figure 
3.8 depict simulation results of the case without disturbances. 
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input u(t) 


0.8 - 
0.6 - 
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Figure 3.7. The Simulation Result of the Control Input. 


reference output and plant output 
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LI Adaptive Control with Bounded Output Disturbance 


Consider a system 


y{t) = K-^u{t) 

p-1 


where 0.5 <C< 3 , and the referenee model is 


p + \ 


where 


Then 


D{p) = p + \. 


B{p) = \ ,A(p) = p-\ ,D{p) = p + \ 


(3.115) 


Assuming = 2 as the true value, the eontroller ean be expressed as 


0.5u{t) = - kg-k^ k^ 


u{t) 

tu 

<t>. 


+ v{t) = ^{t) + v{t) 


(3.116) 


where 


^^=[0.5 2 6f 


= 


1 


1 


u(t) - -u{t) - -y(t) 

p+5 p+5 


(3.117) 


In aeeordanee with Condition 1 and Condition 2, terms are ehosen as 0)^=5, 
(y, = 10 , (Ug = 100, and L = 2> N-M =\. As in the previous example, Figure 3.9 and 
Figure 3.10 show that the stability eonditions are satisfied. 

With appropriate substitutions (p) and Cj,{p) ean be expressed as; 
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IVbgnitude (cB) Phase (cfe^ 


Cp{p) = 


CO,, 


yp + co^) 


CO, 


(p + ®l) 


f 5 ) 

( 7 ^ 

10" 

[p + 5j 

[(p + lof j 


Ce(p) = 


co„ 


{p + COj,) 


100 

(p + 100)' 


bode of (1-(A(s)/(B(s)D(s))) 



Figure 3.9. 


c 

The Bode Plot of 1 

V 


A{p) 

B{p)D{p)^ 


bode of Hyy(s) 



Figure 3.10. The Bode Plot of H^^{p) for 100 Plant Realizations. 
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Figure 3.11 and Figure 3.12 show the results of this simulation. 


input u(t) 



Figure 3.11. The Simulation Result of the Control Input. 



Figure 3.12 The Simulation Result of the Reference Output and Plant Output. 
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3. Example with Sensor Noise 

Consider a plant with transfer funetion 


k^-k ^ 

/v/-v /V, 


0 X / \ 0 

A{p) p + a 


where the parameters are within known bounds 


-2 < a < 0, 0.5 < L < 3 


Let the referenee model be 


1 


1 


Dip) p + l 

First, the LPF, C^ip) is ehosen and is of the form 


Cfip) = 


f \ 

f L ^ 

0)y 

yp + O)^ j 



= Coip)Ciip). 


Again, we ehoose L = 2>N-M=\ . From the bound on 


0 < 7 < = — -1 = ——1 = 5 and the Bode plot of 

k... 0.5 


1 -- 


Aip) 


Bip)Dip) 


Figure 3.13 for a number of random realizations of the plant. 


(3.118) 

(3.119) 

(3.120) 

(3.121) 

, we obtain 
is shown in 
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Magnitude (dB) q IVbgitude (dB) 


10 


abs( 1- ( A(jw)/( B(jw)D(jw) ) ) ) 



e 3.13. Frequency Response of 


Mp) 

B{p)D{p)^ 


for 100 Plant Realizations. 


abs( Hyy(jw) ) 



Figure 3.14. Frequency Response of for 100 Plant Realizations. 
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It can easily be seen that for co>3 rad/see all realizations have magnitude smaller than 
one that required by Condition 1. In this example, (Ug =5 rad/see is ehosen so that 
Condition 1 holds for all realizations, and (U, =10 rad/see so that the gain margin of 
Cp{p) is larger than lOlogjo/^^^ =14J5. This idea ean be applied to the eontroller 
design. 

From Figure 3.14, the magnitudes of all realizations are found to be 

smaller than one required to satisfy Condition 2. Both of these results guarantee the 
stability of H^^{p) in Equation (3.40) for all realizations of the plant. 


The responses of the inputs and output signals for this adaptive eontrol system for 
two different values of cOp are shown in Figure 3.15 and Figure 3.16. To ensure a more 
realistie experiment, bandlimited random observation noise is added. The simulation 
results are shown in Figures 3.17 and Figure 3.18. Notiee that, as predieted by the 
theoretieal arguments in the previous seetion, the peak of the transient response deereases 
as the bandwidth cOp inereases. The relationship between co^, cOp, and the transient 


performanee is evident in Equation (3.87). 


<K 


'' CO, A 


V J 


sup (/){t) 


a<T<t 


(3.122) 


In the next seetion, the upper bounds of parameters and transfer funetions is diseussed. 
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input with wE=10 rad/sec 



output with wE = 10 rad/sec 
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Figure 3.15. Frequency Response to Square Wave for cOj, =10 rad/sec. 
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Figure 3.16 Frequency Response to Square Wave for cOj, =100 rad/sec. 
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input with wE = 10 rad/sec 


output with wE = 10 rad/sec 




Figure 3.17. Frequency Response to Square Wave for = 10 rad/sec with 
Bandlimited Random Observation Noise. 


input with wE = 100 rad/sec output with wE = 100 rad/sec 




Figure 3.18. Frequency Response to Square Wave for cOj, = 100 rad/sec with 
Bandlimited Random Observation Noise. 
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4. Example with Modeling Error 

In this case, the adaptive controller presented in the previous section is used, but 
the plant has unmodeled dynamics as represented by 


y(t) = ka 


1 200 
p-a p^+30p + 200 


u{t) 


(3.123) 


In other words, the system is treated as a first order system, while in reality it has two 
extra poles affecting the high frequency only. It is well known that in this situation the 
adaptive controller ean easily become unstable and sensitive to external disturbances. 
However, the proposed controller seems to retain stability and traeking properties. 
Simulation results are shown in Figure 3.19 and Figure 3.20. 
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Figure 3.19. 


Simulations with Unmodeled Dynamics 
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Figure 3.20. Simulations with Unmodeled Dynamics —^- 
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and Sensor 
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IV. LI ADAPTIVE CONTROL FOR NON-MINIMUM PHASE 

SYSTEMS 

A. MOTIVATION AND BACKGROUND 

The direct Adaptive Control technique presented in the previous chapters 
assumes that the system to be controlled is minimum phase, in the sense that the zeros of 
its transfer function have to be on the left half side of the s-plane. 

It is well known that non-minimum phase zeros present very important 
limitations on the performance of control systems. Most important is the fact that, while 
poles can be moved by state feedback, zeros can be replaced only by pole-zero 
cancellation. Since unstable pole-zero cancellation leads to unstable, uncontrollable, or 
unobservable modes, this operation is clearly not feasible. 

Non-minimum phase systems are very important in some applications and 
therefore cannot be ignored. In particular, systems with flexible appendages most likely 
have zeros on the imaginary axis and are therefore non-minimum phase [89]. 

In this chapter, the adaptive control of this class of systems is discussed. The 
non-minimum phase limitation is avoided by controlling an augmented plant obtained by 
the addition of a system in parallel to the plant to be controlled. This augmenting 
network has to be designed such that its output is zero in steady state, and the overall 
combination is minimum phase. 

Although the control presented for minimum phase systems is restrictive, it does 
exhibit desired tracking in the presence of wide plant uncertainties. Several computer 
simulations show the effectiveness of such a controller. Furthermore, implementation on 
an experimental flexible mechanical arm, with uncertain dynamics, shows that this 
approach can be a viable technique for the control of this class of systems. 


B. FUNDAMENTAL PRINCIPLE AND ANALYSIS OF THE MODIFIED LI 
ADAPTIVE CONTROL 

Along similar lines as in the previous chapters, consider a Single Input Single 
Output (SISO) system with dynamics 
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A{p) 


(4.1) 


with possible non minimum phase zeros. When this is the case, as shown in Figure 4.1, 
we augment it by a parallel system so that the combination system becomes 


Mp) D^(p)) A(p)D^(p) 



D^ip) 


following requirements: 

a. ) the numerator of the combined system (^A{p)E^{p) + kf^B{p)D^{p)^ is Hurwitz; 

b. ) in steady state the output of the augmenting system goes to zero. 

To satisfy Requirement b.) when the input is a constant, Ei(p) is expressed as 

E,{p)^p.E'{p) (4.3) 
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so that the derivative action (the term “ p “) makes the steady state value zero when the 
input is a constant. To satisfy Requirement a.), notice that the numerator of Equation 
(4.2) is the characteristic polynomial of the system shown in Figure 4.2. 


u(t) 


v* / 


A(p) 

-^ 

KE{p) 


i 

L 

E,{p) 


A(p) 


-P* 






yit) 


Figure 4.2. The Other Expression of the New Closed Foop Dynamics with a Designed 

PID. 


Therefore, the augmenting network 


Eiip) 

Diip) 


must be such that its inverse 


Diip) 

E,{p) 


yields a 


stable closed loop system. 


Furthermore, since it is desired that 


Eiip) 

E>i{p) 


is strictly proper. 


E>iip) 

Eiip) 


will not be proper. 


A possible choice is that 


E>iip) 

E,{p) 


is a PID controller, with the integral action to 


satisfy Equation (4.3), and the derivative action to make its inverse proper. For example, 
let a nominal plant be expressed as 


E{p) _ p"+l (P + 17)(P-17) 

A(p) ^ p'^+2p^ ^ p^ (^p+ 1.4142 j)(^p-lAl42j) 


Observe that there are two zeros located on the imaginary axis and therefore, it is non¬ 
minimum phase. 
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In order to make the system minimum phase, a PID eompensator must be 

E,{p) 


designed so that the system in Figure 4.2 is stable. The augmented minimum phase 
system is shown in Equation (4.2). 

From Figure 4.1, the new output is shown to be 


y(0 = y(0 + -^7^w(0 = 

D,{p) 


Mp)Eiip) + B(p)D^ (p) 

Mp)E>iip) 


(4.5) 


Based on the requirements listed above, the degree of E^(p) has to be less thanDj(p). 
Further, according to the final value theorem [89-90], the numerator E^(p) must have 
roots at its “zeros“. Thus define 

E^ip) ^ P'E!{p) (4.6) 

Dj(p) = H- \-ajy (4.7) 

and choose 

E^{p) = p (4.8) 

D^{p) = p^ + p+ 0.\ (4.9) 

Then the transfer function of the augmented system becomes 

BB(p) ^ A{p)E^{p) + B{p)D^{p) ^ p"+pV3p^+l.lp^ + p + 0.1 ^ 

AA{p) A{p)D^{p) p^ + p^+2.1p^+2p^+0.2p^ ^ ' 

with zeros at 


p = [-0.3347 + 1.5021i -0.3347 - 1.50211 -0.1107 + 0.61141 -0.1107 - 0.61141 -0.1094] 


(4.11) 


Which are strictly located in the FHS s-plane. This leads to the augmented system 
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(4.12) 


y(t) = ko 


BB{p) 

AA{p) 


u{t) = 


p^ + p"^+3p^+l.\p^ + p+ 0.1 
/ + p' +2.1p" +2/ +0.2p' “ 


Now, the reference model is chosen as in Equation (3.4): 




1 

P+\ 


v{t) 


(4.13) 


where D{p) = p + \ is a Hurwitz polynomial representing the desired closed loop 
dynamics. The degree of D{p) is N-M , which is equal to the relative degree of the 
augmented plant, and v(t) is an arbitrary external bounded input. 


Designing the standard adaptive controller for the augmented minimum phase 
system follows. In particular, let the filter be chosen as 


Cp{p) = 


f \ 

( 2 ^ 


^ 10 ^ 

f 0 ^ 

10" 

V E + <^o y 



Le + 1oJ 

[(p+iof j 


(4.14) 


It is possible to verify that this guarantees stability of the reference model. 


max max H ( jco) < 1 


(4.15) 


where H^i^p) is as defined in Chapter III. Simulation results shown in Figures 4.3-4.5 
show the stability and asymptotic tracking of the desired trajectory. 
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Figure 4.3. Modified Plant Output and Original Plant Output. 
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Figure 4.4. Simulation Result of Modified LI Adaptive Control System. 
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Figure 4.5. Modified Plant Output Tracks the Reference Output. 


C. SIMULATIONS OF MODIFIED LI ADAPTIVE CONTROL WITH NON¬ 
MINIMUM PHASE SYSTEM 

In order to verify the robustness and stability of the modified Ll adaptive 
controller with PID element, a modified controller is adapted to perform under realistic 
conditions. 


1. Modified Ll Adaptive Control System with Bounded Output 
Disturbance 

Consider a system 


y(0 = K 


Bjp) 

A{p) 


u(t) = 


0.06023p"+0.13378 
/+2.57184p" 


u(t) 


(4.16) 
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where 0.5 <k,<3 , assuming kg=l and the referenee output 


p+\ 


Then 

5(p) = 0.06023p'+0.13378, A{p) = + 2.57mp\ Dis) = p + 1 (4.17) 

This plant has two zeros, one loeated at (0 + 1.49041) and the other at (0 - 1.49041), both 
loeated on the imaginary axis, and therefore the system is non-minimum phase. 

E (p) 

Next, a proper PID eontrol ——- is defined with 

M(p) 


and 


E^{p) = p 


(4.18) 


Dy{p) = 10p^+5p + 0.1 


(4.19) 


It is easy to verify that this eontroller 


E,{p) 

Diip) 


stabilizes the system. 


After the system is 


augmented with 


Eiip) 

E>iip) 


the response is 


-(t)-k - P' +2.8730pVl.3438pV0.6689p +0.0134 ^ 

\ y 0 4 4 / \ \ y 6./x^ 4 . -t r\ ^ r\ 3 . r\ r\^ n ^ 2 \/’ 


AA(p) 


p'’ +0.5p^ +2.5818p" + 1.2859p^ +0.0257p^ 


Now, the referenee model is ehosen as 

y„it) = -^vit) 
p + 5 


(4.20) 


(4.21) 


where D(p) = (p + 5) is a Hurwitz polynomial representing the desired elosed loop 
dynamies. The degree of D{p) is N -M , whieh equals the relative degree of the new 
plant, and v(t) is an arbitrary external bounded input. 
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Furthermore, in order to satisfy the stability conditions of Chapter III, Oq = 30 


rad/sec, <», = 30 rad/sec, and = 50 rad/sec are chosen. Then the LPF becomes 


Cfip) = 


f \ 

( 2 ^ 

nij 


^ 30 ^ 

f 9 A 

30" 

V P + y 



Lp + 3oJ 

l(P+30fJ 


(4.22) 


and 


Ce{p) = 


COr 


{P + (i>eY 


50" 

ip + 50 f 


We can verify that, with this LPF, 


max max H ( ja>) < 1 

a> o<r<r„Ax 


(4.23) 


(4.24) 


and the transfer functions :Hy^(^p), can be obtained 

from Equation (3.42) as before. 

Figures 4.6-4.9 show the results of this simulation. From Figure 4.6, it can be 
seen that the augmented plant is stabilized. In the presence of a bounded output 
disturbance, the settling time is a little longer compared to the example with no 
disturbance in last section. However, due to the effect of the additional PID element, the 
original output tracks the modified plant output and it is still stable, despite the output 
disturbance. In addition, from Figures 4.7 and 4.8, it is obvious that this modified LI 
adaptive control system satisfies the specified requirements of this experiment. It is 
evident that there is no large control action or poor transient performance during the 
duration of execution, and the modified plant output did track the reference output. 
Figure 4.9, shows that the adaptive error is small and bounded, which satisfies the 
requirement of a stable adaptive control system. 


69 



modified plant output and original plant output 

0.5 


0.4 


0.3 

0.2 

0.1 

0 

- 0.1 

- 0.2 

0 100 200 300 400 500 600 

sec 

Figure 4.6. Modified Plant Output and Original Plant Output with Bounded Output 

Disturbance. 



input u(t) 


10 

0 


-10 


-20 



0 100 200 300 400 500 600 


0.3 
0.2 ^ 
0.1 
0 

- 0.1 

0 


sec 

modified plant output 


100 


200 


300 

sec 


400 


500 


600 


Figure 4.7. Simulation Result of Modified LI Adaptive Control System with Bounded 

Output Disturbance. 
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Figure 4.8. Modified Plant Output Tracks the Reference Output with Bounded Output 

Disturbance. 
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Figure 4.9. The Adaptive Error of Modified Plant Output with Bounded Output 

Disturbance. 
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2. Modified LI Adaptive Control System in the Presence of Unmodeled 
Dynamics 

In this case, the same adaptive controller is used as in the previous example, but 
the plant has unmodeled dynamies as represented by 


, 0.06023pV0.13378 

y(t) = K —- 


300 


p"+2.57184p" /+10p + 300 


u{f) 


where A:o=0.1, =0.9* =0.09 <^0 <1.1*^0 =0.11 is assumed. 


(4.25) 


In addition. 


from Equation (4.25), it has two extra poles affecting the high frequencies only. In 
general, in this situation, the adaptive controller can easily beeome unstable and sensitive 
to external disturbanees. Also, the plant has two zeros, one loeated at (0 + 1.4904i) and 
the other at (0 - 1.4904i), both on the imaginary axis and therefore the system is non¬ 
minimum phase. 


Next, the proper control is defined as 


Drip) 


, the inverse of a PID eontrol, where 


Erip) = p 


(4.26) 


and 


Di(p) = 10p^-i-5p-i-0.1. 


(4.27) 


It is easy to verify that this controller 


. , E.(p) , 

system with — - , the response is 


Ei(p) 

Drip) 


stabilizes the system. After augmenting the 


Di(p) 


%,) = Mr) „(,). /+0.6023/ + 2.8730/.1.3438/+0.6689;,.0.0134 


AA(p) 


/ +0.5p" +2.5818p" +1.2859p^ +0.0257p^ 


(4.28) 


Now, the referenee model is ehosen as 
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(4.29) 


with 


yM = -^v{t) 
p + 5 


D{p) = p + 5 (4.30) 

where D{p) = [^p + 5) is a Hurwitz polynomial representing the desired elosed loop 
dynamies. The degree of D{p) is N-M , whieh equals the relative degree of the new 
plant, and v(t) is an arbitrary external bounded input. 


Furthermore, in order to satisfy the stability eonditions of Chapter III, co^ =30 
rad/see, o, = 30 rad/see, and cOj, = 50 rad/see are ehosen. Then the LPF beeomes 


Cfip) = 


f \ 

( 2 ^ 


^ 30 ^ 

f 0 ^ 

30" 

V P + y 

{{p + 0),f j 


Lp + 3oJ 

l(P+30tj 


(4.31) 


and 


Ce{p) = 




{p + (o^) 


50" 

(p + 50)" 


(4.32) 


It is easy to verify that, with this LPF, 


max max H (jco) < I (4.33) 

w o<r<rMAx 


and the transfer funetions :H^^(^p) , H and can be easily 

obtained from Equation (3.42). 

Figures 4.10-4.12 show the simulation results of this experiment. From Figure 
4.10 and Figure 4.11 it is obvious that this modified LI adaptive control system satisfies 
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the specified requirement. It is evident that there is no large control action or poor 
transient performance during the duration of execution, and the modified plant output did 
track the reference output. In Figure 4.12, the adaptive error is shown and is small and 
bounded, which satisfies the requirement of a stable adaptive control system. 
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D. APPLICATION TO A SATELLITEWITH A FLEXIBLE ARM 


In this section, the problem of applying the proposed adaptive controller to the 
control of a flexible arm is addressed. This is part of the equipment available in the 
Satellite Research and Design Center (SRDC) at the Naval Postgraduate School. The 
SRDC aims to recreate realistic space conditions in a laboratory environment. 

The particular equipment used in this experiment is the Flexible Spacecraft 
Simulator (FSS) shown in Figure 4.13. The FSS simulates a flexible arm moving in the 
plane. 



Figure 4.13. Flexible Spacecraft Simulator (FSS). 
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It was developed to experimentally verify the effeetiveness, robustness and stability of 
designed eontrol laws. It is eomprised of a rigid eentral body and a flexible appendage. 
The rigid eentral body represents the main body of the spaeeeraft, and the flexible 
appendage represents a flexible antenna support strueture [91]. In addition, the flexible 
appendage is eomposed of a based beam eantilevered to the main body and tip beam 
rigidly eonneeted to the base beam at a right angle. Five air pads on a granite table 
support the FSS and its flexible appendage to minimize frietion during motion. 

The flexible model is similar to the mass-spring-damping model. Based on 
experimental, real data from FSS, the equations of motion of the flexible model ean be 
developed and are written as 


M {x} + 


(4.34) 


and 


M = 


19.2253 

-1.6190 


-1.6190 

1 


,K = 

“0 

0 

. = 

“O' 

0 

2.2212 


1 


, and X = 


e 


where 6 is the angle of the spaeeeraft hub in radians and q is z. state variable 
representing appendage deformation using modal eoordinates. From Equation (4.34), the 
state spaee dynamies are derived and are written as 


where 


[ Kt) = ApX{t) + B^u{t) 

1 y{t) = CpX{t) 


0 0 10 
0 0 0 1 
0 -2.166 0 0 
0 -2.5718 0 0 
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B 


p 


0 

0 

0.06023 

0.09751 


and 


C,=[l 0 0 0] 

where y{t) = 6{t) is the measured angular position. 


With appropriate substitutions, the transfer function of the plant becomes 


B{p) 0.06023p"+0.13378 

3^(0 = =- 4 ^ _ 7 ^(0 


A{p) 


pV2.57184p' 


(4.35) 


It can be seen that, this flexible model has two non-minimum phase zeros located on the 
imaginary axis. The root locus plot of this plant is shown in Figure 4.14. 


Root Locus 



Figure 4.14. The Root Locus Plot of the Non-minimum System 

0.06023p"+0.13378 
y{t) = —s^— u{t) . 

/+2.57184p" 
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It can also be verified that the PID eontroller with transfer funetion 


D^{p) _ 10 /)^+ 5 /) + 0.1 
E,{p) p 


(4.36) 


stabilizes the system. Also, it provides suffieient margin so that parameter perturbations 
do not affect the stability. Then, the nominal augmented system is 


y{t) = 


B{p) , E^{p) 


- + 


A{p) 


u{t) 


(4.37) 


and beeomes 


_ - BB{p) _ + 0.6023/+2.8730pVl.3438pV0.6689p + 0.0134 ^ 

^ y 0 t t / \ ft . r\ n 5 . r\ n (~i -t (~i 4 . -t r\ n n /'\ 2 


AA(p) 


^”+0.5^+2.5818^ +1.285V +0.0257p^ 


(4.38) 


where kg=l and A:o = 0.1. 


Sinee the relative degree is A - M = 1, the referenee model is ehosen as 

p + 5 


(4.39) 


with v(t) representing an arbitrary external input. It ean be easily verified that the 
augmented system is minimum phase, and ean be adaptively eontrolled to follow the 
trajeetory of the reference model of Equation (4.39). 

Proeeeding as in the previous ehapter, we ehoose the fdters Cf{s) and Ci,{s) for 
the purposes of eontrol and adaptation. With co^ = 30 rad/see, co^ = 30 rad/see, and 
CO=50 rad/see ehosen, the EPF beeomes 


Cfip) = 


c \ 

f 2 ^ 

CO^ 


r 30 ^ 

f 9 A 

30' 

V P + y 

j 


Ip + 3oJ 

l(p+3o)y 


(4.40) 


and 
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L 


(4.41) 


Ce(p) = 


a>^ 


{p + (i>E) 


50 ^ 


From the plot of in Figure 4.15, it can be see that the magnitude is 

always smaller than “one”, thus guaranteeing stability of the transfer functions ; (p), 

^yo{p)^ ^„o(;^)>and in Equation (3.42). 


abs(Hyy(jw)) 



Figure 4.15. The Plot of the |//^y( 7 (»)| 


Figures 4.16 and 4.17 show the results of this experiment. Figure 4.16 shows the 
trajectories of the reference model, original model, and the modified plant. Figure 4.17 
shows the results of control input, reference output, modified plant output, and original 
plant output, respectively. From these figures, it can be seen that the modified plant with 
PID control results in the plant becoming stable. Comparing the figure of the reference 
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output, it is evident that the modified output after adding the PID does track the trajectory 
of the reference model. In addition, no large control action or poor transient performance 
is exhibited during the entire execution period. This result validates that the modified 
adaptive controller with PID design does solve the flexible problems in non-minimum 
phase systems. This is the reason that such a controller is implemented in this 
experimental adaptive system. 



Figure 4.16. Trajectories of the Reference Model, Original Plant, and Modified Plant 
with Modified LI Adaptive Controller and FSS. 
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Figure 4.17. Experimental Results with Modified LI Adaptive Controller and ESS. 
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V. CONCLUSIONS 


A. SUMMARY AND CONTRIBUTIONS 

In this dissertation, an Adaptive Control system has been redesigned for Single 
Input Single Output (SISO) systems to improve both transient response and robustness to 
external disturbanees and unmodeled dynamies. In partieular, it was shown how two 
properly designed filters, one in the eontrol loop and one in the parameter estimation bloek, 
yield a transient response due to adaptation with an arbitrarily small maximum norm. Also, 
based on a number of simulations, it was verified that the system remained stable in the 
presenee of external disturbanees and unmodeled dynamies. 

This approaeh was extended to non-minimum phase systems by assuming some 
stability eonditions using a non-proper eompensator, sueh as a PID eontroller. In the latter 
ease, the effeetiveness of this approaeh was demonstrated by applieation to a system with 
flexibilities, having a non-minimum phase transfer funetion. 

The eontributions in this dissertation are as follows: 

1. For minimum phase, SISO systems global stability has been established 
together with eriteria on the ehoiee of the various Alter parameters. Also, 
upper bounds of the transient response due to adaptation have been 
eomputed and shown how they are affeeted by the filters’ bandwidths; 

2. For non-minimum phase systems, we have shown that the use of an 
augmenting system in parallel with the plant makes the overall system 
minimum phase, so that the Adaptive Controller can be applied. The 
criterion of choosing the augmenting network is by the inverse of a PID 
controller, which stabilizes the plant. Although the PID controller has a non¬ 
proper transfer function due to the derivative action, its inverse is proper 
thus easily implementable in the presence of additive disturbances. The 
augmenting network is chosen so that it has no effect in steady state, only 
during the transient response of the system; 

3. Application to a very lightly damped flexible system with two degrees of 
freedom on a plane shows the effectiveness of this approach in an actual 
implementation. 
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B. FURTHER WORK 


The study in this dissertation presents a preliminary investigation on an Adaptive 
Control system with the intent of making it attractive to engineering problems. There are 
still a number of issues that need to be addressed, namely: 

1. Robustness in the presence of external disturbances and unmodeled 
dynamics. Although simulation results show that this controller still keeps 
its stability properties with modeling errors, a need to establish how mush 
uncertainties can be tolerated while still guarantee stability has to be 
addressed. In particular the claim of arbitrarily small transient errors due to 
adaptation has to be reformulated to take modeling errors into account. The 
expected final result will be a usual tradeoff between performance and 
robustness; 

2. The extension of this algorithm to non-minimum phase systems seems to be 
very important in view of the fact that systems with lightly damped 
flexibilities are very hard to control. The algorithm presented assumes that a 
stabilizing non-proper controller (such the PID) is known to the designer. 
The issue to be further investigated is the best way of determining this 
stabilizing controller, in view of the given uncertainties of the system to be 
controlled. From an engineering standpoint this would correspond to the 
combination of two controllers: a standard fixed PID control just to stabilize 
the system and the adaptive controller to improve its dynamic performance. 
These two controllers (fixed and adaptive) interact with each other, but their 
effects on the system performance and robustness need to be further 
investigated; 

3. Finally most systems of interest have Multiple Inputs and Multiple Outputs 
(MIMO). This is the case of the flexible arm presented in the last section: 
extension to an actual implementation in space would require dynamics in 
all three dimensions with actuators in all three axes. For this sort of 
problems there are two possible ways of extending the theory and design. 
One way is to take the whole dynamics of the system into account and 
extend the SISO approach to MIMO. Along these lines a considerable 
amount of work has been done for standard Adaptive Control systems and it 
can be extended to the LI approach. However in some cases as in the 
flexible arm in space, the MIMO system can be modeled as a set of SISO 
systems with cross coupling dynamics. In this case it might be easier to 
design separate independent SISO controllers and treat the cross couplings 
as disturbances. In this way, if the SISO controllers have sufficient degree 
of robustness the overall system will still be stable. These show important 
tradeoffs of this class of design. 
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APPENDIX A. 


Equation (3.31) can be written as 




X ,(0 


cr 


-2 L- 


VV,(0| +\\X{t)\ 




X ,(0 


cr w^{t) + X{t) 


-e{t) 


where ju - ju. (t) > 0. 


(A.1) 


Define the positive define funetion 


v{e{t)) = U{tye{t) (a.2) 

then 

vm)) = e\tym (A.3) 


and 


vm))=-M- 


e^{t)x{t) 


2N 


w^it)^+ X(t) 


-«’(0 + X ( a - a ,( 0 )- 


xyt) 


cr 


-2 L- 


VV,(0| +||A(0|| 


-e{t). (A.4) 


The rightmost term is nonzero only when 


0, {t) = e- ^ < 0 and A,. {t)e{t) > 0 


(A.5) 




or 


L ^ (0 = ^ - ^ 0 and A,. {t)e{t) < 0 

whieh makes it always non-positive. 

Substituting for the error term e{t) = 0^ {t)X{t) + Wj:{t) the following expression is 
obtained 




-A 


^(0 + ^? 7(0 




(A.6) 
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where 


r 


s{t) = 


e^{t)x{t) 


<j^ WgCO + x{t) 


< 


77(0 = 


Wsit) 


.s{t) = ■ 


e^{t)x{t) 


V. 


^a^\w^{t)f +\\X{t)f +\\X{t)\ 


(A.7) 


Now, with any t, T > 0 and integrating Equation (A.6) between t and t+T 


t+T 




t+T 

dT+— J |;7(r)|^ Jr (A.8) 


yields 


t+T 




dr < 


2// 4 


t+T 


j ri{T)^ dr. 


(A.9) 


Definition of rj{t) in Equation (A.7) yields the upper bound on the rightmost term as 


t+T t+T 

j dr - 


0 if W£(r) = 0 for all r 


(A.10) 


vvglr) j 

dT< ^ 


r o-^ w^{t) 


otherwise. 


cr 


The final results beeome 

a. Equation (3.33a), eome easily sinee the initial estimates are within the 
bounds 

^ ^(0 ^ 0i.MAX ’ ^nd 


^, (0 ^ 0i,MAX ^ ^ 0 

4(o<^,„^4(o>o 


(A. 11) 
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b. From the definition of £{t) in Equation (A.7) 


6*' 


1 

(0X(0|<|^(0|(^^|w.(0r+||^(0ff 


(A. 12) 


where we ean bound 


j, 


1 r 


£(^) + -ri(T) 


dr 


+ 


i : 


\Hr) 


1 

2 ^2 

Jr 

y 


(A.13) 


Finally, eombining Equations (A.9), (A.12), and (A.13) yields Equation (3.33b) 
with bound Equation (3.34). 

e. Sinee 


\\X{t)f\e{tf 


0{t) 


<A 


<A 


e{t)' 


(a>,(0r+||A(0f) +||A(0 


(A. 14) 


and substituting for e{t) = 6^{t)X(t) + Wj,{t) yields 


e{t) 




( 1^1 

= £{t) + ri{t)= £{t) + -ri{t) +-ri{t) 
V 1)1 


(A.15) 


This expression has the same upper bound as Equation (A.13) so that Equation (3.33o) 
follows easily with bound again of Equation (3.34) 


QED 
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APPENDIX B. 


Proof: From Equation (3.20), the expression for the plant output y{t) is 


y{t) = 


-kju{t) + + k'^^yit)) + w{t). 

Dip) 


(B.l) 


Combining Equations (3.21), (3.22), and (3.23), and rearranging Equation (3.20) as 
follows yields 


0^^{t) = Dip) (y(0 - wit)) + (^0 - kjuit) - kf^uit). (B.2) 


Erom Equation (3.1) 


K^it) = ^^{yit)-wit)) 


(B.3) 


Substituting Equation (B.3) into Equation (B.2) results in 


0^<Pit) = Dip)yit) - Dip)wit) + iko - kjuit) - yit) + it). 

Bip) Bip) 


(B.4) 


hen 




B^^it) = Dip) 


Aip) 




V Bip)Dip) 


yit) + ikg - k^)uit) - Dip)wit) + it). 

Bip) 


(B.5) 


Substituting Equation (B.5) into Equation (B.2) yields 




r 


Dip) 
(1-C,(p)) 
Dip) 


Aip) 




V Bip)Dip) 


yit) 


ikf, - kjuit)-{\-Cpip)) wit) + Woit) + wit) 

Dip) Aip) 


(B.6) 
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kji{t) = Cp {p)[e'^ + v(0) - Cp{p)D{p) 


1 -- 


A{p) 


B{p)D{p) 


-Cp {p){k^ - kju{t) + Cp {p)D{p)w{t) - Cp (p) 


Mp) 

B(p) 


y(t) 


Woit). 


(B.7) 


Equation (B.6) can be adjusted as 

j(0 = (i-c^(e)) 


\ Mp) ^ 

V B{p)D{p)j 


y( 0 ++ XT +v( 0 ) 


CFiP)rnT. 


-w{t) + Cp {p)w{t) + [\-Cp{p)) 


D(p) 

^ Mp) ^ 

B{p)D{p) 


D(p) 


w^{t) + w{t) 


then 


y{t) = {\-Cp{p)) 


\ Mp) ^ 

V B{p)D{p)j 


y(0 +(0^(0+v(0) 


CFiP)rnT. 


+ {^-Cp{p)) 


( 


A{p) 




Bip)D(p) 


D(p) 


wJt) + Cp{p)w{t) 


D(p) 


(B.8) 


and finally 


k^ (1 + rCp{p))u{t) = -Cp{p)D{p) 


1 -- 


A{p) 


B{p)D{p)^ 


-Cp(p) 


Mp) 

B{p) 


wM)-B>{p)w{t) 


y(t) + Cp (p)(M (t)<p(t) + v(0) 

(B.9) 


Then k^u{t) is isolated; 


k^u{t) = 


-Cp{p)D{p) 


{\ + yCp{t)) 

CpiP) 


\ A(p) ^ 
V B(p)D(p) 




(l + yc^(0) 


Mp) 

B{p) 


wM)-B>{p)w{t) 


(B.IO) 


Substituting kji{t) from Equation (B.IO) into Equation (B.8), results in the proof of 
Equations (3.41) and (3.42). 
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where 


y{t) = Hyy{p)y{t) + H^o{p)CFiP)[^^ (0^(0 + v{t)) + H ^^[p)w^ (t) 
k^u{t) = H^y{p) y{t) + H^^{p)Cp {p)[0^ it)(p{t) + v(0) + H^„{p) w^{t) 


(3.41) 


and y 


H„(p)=(i-cyp)) 


A(s) 


B(p)D(p) 


1 + 


y 


^vo(i^) 


H^P) 


(1 + /) 

(i+rCp{p))D(p) 

(l-Cy(p)) f A{p) ^ k{p) ^ 
(l + yCp{p)){B{p)D{p) P{p)D{p) ^ 


hM 

Huy{p) 

hM 


1 

\+rCp{p) 

CF{p)D(p) f A{p) ^ 

\ + yCp{p) yB{p)D{p) ^ 

Cf{p)D{p) ( Ajp) ^ kjp) 
i+rCp{p) {b(p)D(p) p(p)D(p) ^ 


k{p) 

P{p)D{p) 


k k 

— -1 in the interval 0 < / < y,^^^ = —^ -1. 

K K 


(3.42) 


QED 
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APPENDIX C. 


Let 

a. C^{p) be the transfer funetion of an exponentially stable, eausal Linear Time 
Invariant (LTI) system, 

b. d(t) be a bounded veetor 

e. (p(t) e ,e{t) e i? be the input and output of the following Linear Time Varying 
operator 




(C.l) 


then 


t t 

e(t) = I C^(t -T)d^ {T)^(T)dT C^(t -T)d^(t)^{T)dT (C.2) 

0 0 



Figure C.L Region of Integration for Equation (C.2). 
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Since 


e{t) = e{Qi) + \^{X)dX (C.3) 

0 


Substituting Equation (C.3) into Equation (C.2) yields 

t t 

e (?) = - J J Cg (? - r)(9(/l)^zi(r) J/lJr (C .4) 

0 T 


It is easy to see from Figure C. 1 that the region of integration can be expressed as 


0<X<t 

0<T<X 


(C.5) 


which yields 

' . ^ 

e{t) = -J 0^(X)j c^(t-T)^(T)dTdX (C.6) 

0 0 


A simple change of integration variables yields the result. 


QED 
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